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Goals 


Assimilate satellite retrievals of soil moisture into a regional (3-km) 

land surface model (SPoRT LIS running Noah 3.3). 

- Take advantage of high-resolution geophysical properties, best 
available atmospheric forcing, and latest satellite measurements 
on soil moisture 


Predicted impact 
- Improved representation of fine-scale soil moisture fields 


- Better depiction of gradients and structure for coupling with 
NWP models at convection-allowing resolution (3 km) for 
regional weatherforecasting 

Demonstrate impact on: 

— LSM soil moisture field 
— coupled NWP forecasts 
Transition a real-time version of LIS output to end users. 


Use experience from SMOS assimilation to implement SMAP. 
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Motivation 


— Improve model depiction of soil moisture and related variables 


(Direct Applications) 
drought monitoring, flood forecasting, agriculture 





— Better numerical weather forecasts using coupled NWP/LSM 
Available moisture affects humidity, sensible/latent heating, diurnal heating rate, 
and convection. 
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‘Short-term Prediction Research and 
Transition (SPoRT) Center 


Mission: Transition unique NASA and NOAA observations and 
research capabilities to the operational weather community to improve 
short-term weather forecasts on a regional and local scale. 

"= Close collaboration with numerous WFOs and 

National Centers across the country 

= SPORT activities began in 2002, first products 
to AWIPS in 2003 

= Co-funded by NOAA since 2009 through 

Proving Ground activities 


= Proven paradigm for transition of research 
and experimental data to operations 
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Benefit: 
= Demonstrate capability of NASA and NOAA (nt tpt Pa 


© National C ester I -vahustion Partner 


experimental products to weather applications  ¢ wocs--<rm» 
and societal benefit 


" Take satellite instruments with climate 
missions and apply data to solve shorter-term 
weather problems 
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Land Information System (LIS) 












Land Surface Models (LSM) 
(15 min-1 hour) 


Surface Fluxes: 


Topography, 


Soils 
(Static) 


Onet,Qle.Qh,Qg 
Swnet, LWnet 





Noah, CLM, VIC 









Soil States: 






Vegetation, 







LAI Moisture, 
tL 
Temperature 
Monthly 
2 Profiles 














Water Balance: 
Runoff, Drainage, 
Baseflow (some) 


Meteorological 
“Forcings” 
(Hourly-3 hourly) 









Surface States: 
Snow 
Skin Temperature 
LAI (some) 










References: 
Kumar et al. (2006) 
Peters-Lidard et al. (2007) 








and Surface States 
(Variable) 






Framework for running LSMs incorporating a wide variety of 
meteorological forcing data and land surface parameters 
=" Developed by NASA-GSFC 
= Includes data assimilation capability. 
=" Can be run coupled with Advanced Research WRF. 
Using Noah 3.3 Land Surface Model (LSM) within LIS 
SPORT maintains near-real-time and experimental LIS runs 
= SE US (3-km), shared with WFO’s 
" East Africa, shared with Kenya Meteorological Service (KMS) 
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SPoRT-LIS total column soil 


moisture displayed in AWIPS Il 
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East Africa LIS domain 





For more on SPORT LIS, see also: 


Real-time Land Information System over the Continental U.S. 
for Situational Awareness and Local Numerical Weather 
Prediction Applications 
(Case et al., Hydro 3.3) 


Full Continental U.S. (CONUS) domain with 
0.03° (lat/lon) grid resolution 


Unique characteristics of SPoRT-LIS: 


— Real-time S-NPP/VIIRS Green Vegetation Fraction 
— Albedo scaled to input vegetation 


— Restart simulation strategy to produce real-time 
output (timeline below) 
SPORT-LIS ingested and displayed in AWIPS II 
at select NOAA/NWS weather forecast offices 


Land surface variables available to initialize 
modeling applications (WRF and STRC/EMS/UEMS) 





Current SPoRT-LIS CONUS domain, 
as displayed in AWIPS Il 


SMOS and SMAP 


¢ L-band radiometers (and radars) 
can be used to estimate soll 
moisture near the surface 


— Compared to higher frequency 




















instruments: 
o Sees deeper in the soil (“1-5 
cm) 
o Better vegetation penetration [SSS aS ES 
o Higher sensitivity (accuracy) Agency NASA/ ESA NASA 
. : JAXA 
¢ SMAP radar gives improved 
; ‘ Launch | 2002 | 2009 Jan. 2015 
horizontal resolution 
oo. Orbit | Polar | Polar Polar 
° Assimilating retrievals from Soil Sensor Passive Passive Passive Active Combined 
Moisture and Ocean Salinity Type (Failed limited 
° ration 
(SMOS) satellite july 
2015) 
° Prepa FINE for assimilation OT NASA frequency | 6.9GHz | 1.4GHz 141GHz 12 
Soil Moisture Active/Passive (C-band) | (L-band) GHz 
(SMAP) retrievals Resolution | 56 km | 35-50km 36km 3km 9km 
— SMAP has higher resolution product Accuracy 6cm?/cm? | 4cm3/em? 4cm3/ 6 cm?/ 4 cm3/ 
but due to failure of radar, time period cm® cm° cm? 








is limited to a few months. 


Data Assimilation in LIS 


A One-Dimensional Ensemble Kalman Filter: 
Assimilating an Observation 


Posterior L/ Obs. Likelihood 
Ely. , rt | dor 
Background Parl eb i. 
. 2 4 


Figure from J. Anderson, NCAR. 





e Uses Ensemble Kalman Filter in LIS 


¢ Combines Background (Model) and Observations (Satellite 
Retrievals), weighted by their uncertainties, to provided a 
new analysis 

¢ Observation operator relates the top model layer of soil 
moisture (0-10 cm) to the bias-corrected observations (~5 
cm). 

¢ Better depiction of top layer can improve deeper layers 
through infiltration and diffusion. 


Bias Correction 
¢ Data assimilation systems generally assume unbiased observations. 


¢ In general, SMOS observations (retrievals) are drier than the model but have a 
higher dynamic range. 


¢ CDF-matching is commonly used in land surface modeling (forcing observations to 
match model distribution) 


Innovations (Obs-Bkgd) 























U ted Bias Correction Innovations 
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Bias Correction 


Correction Curves 
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¢ LIS can apply point-by-point correction curves. Many implementations generate 
climatologies of model and obs at each grid point. 


¢ We tested three variations of CDF matching, aggregating spatially to increase sample 
size. 
— Single uniform correction 
— Soil-type based 
— Vegetation-based 
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SMOS Experiment 


Precursor to SMAP Land Cover es 


Southeastern/Central USA 3-km domain ae 


40N +—= 


Mixed Forest 
Woody Savanna 


EB Forest 


MODIS/IGBP Vegetation Type 
STATSGO Soil Type 


Cropland 





. Crops/Nat. Veg 
Daily MODIS GVF bees, 
North American Land Data Assimilation 2 ,,| cee 
(NLDAS-2) forcing nf eevee 





Precip: Stage IV (radar+gauge) 


1-yr spinup, 1 month perturbations, 32 ensemble 
members 
Experiment run March-October 2011 
Control (Open loop with perturbations) 
DA run (3 different bias corrections + no correction) 
Validation 

¢ North American Soil Moisture Database 

¢ Due to scale mismatch, expect correlations to be 

most useful metric 


3-km results (14 May 2011) 
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ANALYSIS 
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OBSERVATIONS 


Increment 12214MAY2011 


44N 
42N 
40N 
38N 
36N 
J4N 
32N 
SON 
25N 
26N 


24N 


ggW SBW gaw gow arw Baw aw 7aw 75W 72H bow 
Soil Moisture (cm?/cm*) 


INCREMENT 





gow ob gS Saw aw BAW aiW FB FO F24 6oW 





0.5 


0.45 


G4 


0.35 


0.3 


Q,25 


G.2 


0.25 


0.05 


Q,025 


—0.025 


—0.05 


—0.1 


—0.16 


—0.2 


—0.25 


5em Soil Moisture 


5em Soil Moisture 


Sém Soil Moisture 


SMOS DA Validation 


BOLUS eU Zz 


1 















































ve pias ErrSD Corr 
FL —0.042 0.134 Fo 
Oo oe 
DA 0.037 0.104 0.917 1.Prairie SP, MO 
0.4 
a ~Anlhaewonnll Monn OUN, IPS, _f 
- NA apr, A= 
ao ANS JH 
0.1 
0.0 ! l l 
Mar 1 Apr 1 May 1 Jun 1 Jul 1 Aug 1 Sep: 1 
Date (2011) 
QO070800212 128 
ae Bide ErrSD Corr | | | = 
0.5E OFL O,055 0,071 0.678 = 
= DA 0.036 0.056 0.886 ° = 
— = 
= 2.Stillwater,OK 4 
03 RAS \ \ N\ AS 
ETA Rn AN YE. oat\Ne ISS IN 
0.2 av varal Pawel OP \A AP ty = 
Oa = 
= = 
= = 
O.OF | | | | ! | = 
Mar 1 Apr 1 May 1 Jun 1 Jul 1 Aug 1 Sep 1 
Date (2011) 
OLOOOS 22 1O/ 
a es ErrSD Corr | | S| 
0.5 OFL TS 0.101 0.702 a 
DA =O. 7S 0.070 0927 : . a 
3. Pittsfield _ 
e ) = 
Se ea ag 
WS ai 
0.1 2 
0.0 | | | | | = 
Mar 1 Apr 1 May 1 Jun 1 Jul 1 Aug 1 Sep i 


Date (2011) 








¢ 0-10 cm model soil moisture 
¢ Compared open loop run to 
SMOS DA run. 


Results from validation against 
soil moisture networks in US 
(North American Soil Moisture 
Database) 

¢ Better correlations 

¢ Improved dynamic range 


SMOS DA Validation 


O— 10 cm, Change In Correlation Root Zone Change in Correlation 
= ——— : a = 75 —90 = EO 75 
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SMOS DA Validation 
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Summary of Bias Correction Results 


Variable 0-10 cm Soil Moisture 

# Stations 194 

Experiment OPL NOBC BC1 BCS BCV 
Bias -0.000£0.011  -0.026£0.011  -0.023+0.011  -0.005+0.011 -0.025+0.011 
RMSE 0.082 0.005  0.087+0.006  0.086-£0.005  0.082£0.005 0.087 0.006 
aa 0.046+£0.003  0.043£0.002  0.043+0.002  0.044+£0.003 0.043 0.002 
Correlation —0.45140.023 —(0.573-£0.027 —-0.569£0.026 (0.539 +0.025 0.561 + 0.026 
Variable Root Zone Soil Moisture 

# Stations 137 

Experiment OPL NOBC BC1 BCS BCV 
Bias 0.038+0.015  -0.013+0.016 -0.002+£0.016 0.014+£0.016 0.009 +0.017 
RMSE 0.093£0.008  0.094+£0.008  0.092-£0.008  0.092-£0.008 0.094£0.008 
a 0.037+0.003  0.040+0.003 0.036£0.002 0.038 +0.003 —:0.038 + 0.003 
Correlation 0.6724 0.040 —0.685-£0.043 0.680 0.043 0.6674 0.042 —-0.677 4 0.045 


Experimental error statistics with 95% confidence intervals for 0-10 cm layer soil moisture, 
verified against Texas A&M North American Soil Moisture Database in situ observations from | 
April to I October 2011. OPL: Open Loop; NOBC: Data Assimilation Only, BCI: single bias 
correction; BCS: soil-based bias correction; BCV: vegetation-based correction. The best 
statistics in each category are in bold font. 


¢ All DA runs improved correlation signficantly in upper zone (0-10 cm). 

¢ Soil type correction did best job of reducing bias (as compared to stations) 

¢ Representativeness error could be reduced in future by comparing against 
COSMIC probes. 


SMAP soil moisture assimilation 
¢ Original plan: assimilate combined active/passive (L2) 
retrievals (9 km) 
¢ SMAP radar failed July 2015 
¢ New plan: assimilate passive (L2) retrievals (36 km) 


¢ Alternative: possible higher resolution products from SMAP 
science team? 


Name SPAN 











Launch | Jan. 2015 

Orbit Polar 

Sensor Type Passive Active Combined 
(Failed (limited 
July 2015) duration) 

Frequency 1.41 GHz 1.2 GHz 

Resolution | 36 km 3 km 9km 





Accuracy 4 cm?/cm? 6cm3/cm? 4cm?/cm? 





Sampling Strategy 


¢ Level 2 data are available on 36-km EASE grid 

¢ To take advantage of high resolution geophysical properties (topography, 
vegetation, soils), running model at 3-km 

¢ SMAP observations are assimilated at each model grid point in their FOV 


LIS grid (3-km) 

















Some QC applied on LIS grid impo ea 
Depends on LSM/variable \ ae ae 
(e.g. Noah3.3+soil moisture) Pics SMAP (passive) 
¢ Precip (changed to had aa 36-km cell 
Rie 
1 mm/hr) pss 
¢ Frozen ground ABBE Data flag-based QC applied 
¢ Snow on ground Bm f at observation resolution 
Sir , 
¢ GVF>0.7 nnn, ¢ Retrieval Quality Flag 
¢ Extreme values (new in PE ¢ Vegetation Opacity 
LIS 7) ¢ Vegetation Water 
e “Forest” land class ¢ Frozen Ground Fraction 


SMAP and LIS grids are not 
aligned. Near boundaries, keep 
only one observation per cell 
(closest good ob) 


Bias correction will be 
applied on LIS grid. 


First assimilation results 


SMAP Soil Moisture Observation SMAP Soil Moisture Innovation (Ob-Bk) 
18Z 1 May 2015 18Z 1 May 2015 
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WRF impact tests (planned) 


LIS/WRF Coupling ¢ Coupled LIS/WRF runs 
—NWP provides forcing for LSM 
Noah LSM in LIS —LSM provides fluxes and surface 





conditions to NWP model 
e Assess impact of SMAP DA on NWP 
for coupled runs 


LIS Output Tome eel —Verify NWP forecasts against surface 
pe lg a aelealy obs, soundings, and precipitation 
analyses 


—Examine impact on significant events 





WRF NWP Model 





Validation Datasets 
Domain | T,q,wind _| Precipitation 
conus. [| MADIS” | MRMS 
‘EastAfrica | WMO network |GPMIMERG 


GFS Atmospheric 
Boundary 
Conditions 








Summary and Plans 


Implemented SMOS data assimilation in Noah LSM within LIS 


¢ Improved correlations with ground observations for upper layer (0-10 cm) and root 
zone (10-100 cm). 


Currently testing SMAP assimilation (passive 36 km L2 product) 


Future Plans 


¢ Validation against NASMD including COSMIC probes (reduced representativeness 
error) using LIS Validation Toolkit 


¢ Implement SMOS/SMAP DA in near-real-time SPoRT LIS product 
¢ Validation of bias correction methodology 
¢ Coupled LIS-WRF experiments using NU-WRF 

— NWP validation over US and East Africa 


— Expect more dramatic improvement over Africa where observing networks are 
less extensive. 


¢ Implement DA in near-real-time SPoRT LIS runs 


— Transition products to NWS and international partners 
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Questions and Comments? 


Contact information: 
clay.blankenship@nasa.gov 
http://weather.msfc.nasa.gov/sport 
Facebook: NASA.SPoRT 
Twitter: @NASA_ SPoRT 


